Correlation functions of the one-dimensional attractive Bose gas 
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The zero-temperature correlation functions of the one- dimensional attractive Bose gas with delta- 
function interaction are calculated analytically for any value of the interaction parameter and number 
of particles, directly from the integrability of the model. We point out a number of interesting 
features, including zero recoil energy for large number of particles, analogous to a Mossbauer effect. 



Recent experiments on trapped one-dimensional (ID) 
atomic gases [1, 2] provide a unique opportunity to study 
the effects of quantum correlations and fluctuations for 
one of the paradigms of strongly-correlated systems: the 
Lieb-Liniger interacting Bose gas [3] . One of the striking 
features of these experiments is that the effective ID cou- 
pling c can be tuned to essentially any (positive or nega- 
tive) value. On the theoretical side the model is exactly 
solvable, and therefore physics beyond mean-field can be 
reliably investigated. Although the repulsive regime is 
well-understood, the attractive case is less commonly dis- 
cussed in the literature. This case is unusual because the 
atoms can bind together: in fact, the ground state of N 
particles of mass m in infinite space is a clump whose 
binding energy is E a = -mc 2 N(N 2 - l)/6ft 2 [4]. This 
clump behaves like a particle of mass given by Nm. 

Our purpose here is to provide a detailed analysis of the 
dynamics of this system, by calculating observable corre- 
lation functions analytically. We give expressions for the 
one-body function and dynamical structure factor (DSF), 
which are experimentally accessible using ballistic expan- 
sion and Bragg spectroscopy, and highlight their impor- 
tant features. The attractive regime has found renewed 
experimental interest with the observation of bright soli- 
tons in quasi ID harmonically trapped condensates [2]. 

The Hamiltonian of the Lieb-Liniger model is given by 
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with c = — c > the interaction parameter, m the mass 
of the particles (atoms), and the sum runs over all pair- 
ings. In terms of experimental parameters the ID cou- 
pling constant is c = h 2 /many, where a\D is the effective 
ID scattering length that can be tuned via Feshbach res- 
onance or transverse confinement [5]. For definitcncss, 
we consider a system of length L with periodic boundary 
conditions. From now on we fix h = 2m = 1. 

Hamiltonian (1) is diagonalized by the Bethe Ansatz 
[3]: the eigenfunctions are superpositions of plane waves, 



A' 



over all permutations P of the 



momenta (rapidities) A. The Ap coefficients are func- 
tions of the two-particle scattering phase shifts obtained 
from (1), and periodicity of the wavefunction requires the 



set of rapidities {A} to be solution to the Bethe equations 
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For the repulsive case, all solutions are real [6], and the 
wavefunctions are scattering states of N atoms. For the 
attractive case c > 0, however, the physics is completely 
different. Complex solutions to the Bethe equations are 
allowed: atoms can therefore bind and form new types 
of particles which are stable under scattering. A general 
eigenstate is therefore made by partitioning the N atoms 
into a set of N s bound-states. The rapidities associated 
to a bound-state of j atoms form a regular pattern in the 
complex plane which is called a string [7]: 



AS a =K + -(j + l-2a) 
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Here, a = 1, j labels the rapidities within the string, 
while a = l,...,Nj labels the set of Nj strings of length 
j. Note that N = YljjNj, with the total number of 
strings given by N s = J^j Nj . 5^ a are deviations which 
fall off exponentially with system size L. Perfect strings 
(i.e. with <5 = 0) are then exact eigenstates in the limit 
L — > oo for arbitrary N. Thus we consider the limit 
L -> co at fixed N [18]. 

Such an eigenstate will have momentum and energy 
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where a ) represents the sum over strings in the eigen- 
state. The set of N s rapidities A^ obey a set of reduced 
Bethe equations obtained by substituting (3) into (2). 

The ground-state corresponds to a single TV-string cen- 
tered on zero. Excited states can be classified in terms of 
their string content. One-particle (or one string N g = 1) 
states are obtained by giving a finite momentum K to 
the ground state. Their dispersion relation is co = K 2 /N, 
meaning that in the large N limit they form a flat disper- 
sion band degenerate with the ground state. Two-particle 
states (N s = 2) are composed of two strings made of M 
and N — M particles respectively. Multiparticle states 
with N s > 2 are similarly constructed. 
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Form factors and correlation functions. We consider 
the zero-temperature density-density correlation function 



SP{x, t) = (p{x, t)p(0, 0)) with p{x, t) = S(x - Xj ), 

whose Fourier transform is known as the dynamical struc- 
ture factor (DSF), and the one-body dynamical cor- 
relation function of the canonical Bose field ^(x,t), 
S*(x,t) = (*t(a;,i)tf(0,0)). By Fourier transform 
O(k,oj) = jf dx f^dte^'-^Ofat), these can be 
written as a sum over intermediate states fi, 



S°(k,uj) = 2ttLJ2 
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where GS denotes the ground state, and the form factor 
(FF) = (n\O(p,0)\GS) = ±(p\0 K JGS) depends on 
the operator O and on the state p {Ok is the Fourier 
transform of the operator O at momentum K). \\p\\ de- 
notes the norm of state p. For the Lieb-Liniger model, 
FFs and norms were calculated for general Bethe states 
in Refs. [9-12]. They are given by the determinant of 
matrices whose entries are rational functions of the ra- 
pidities of the two eigenstates involved (we do not report 
these rather complicated expressions here). 

In the repulsive case, the Bethe equations cannot be 
solved analytically. It is however possible to solve the 
Bethe equations numerically at finite N and L, determine 
the corresponding FFs, and perform the sum over inter- 
mediate states by selectively scanning the Fock space, 
thereby obtaining very precise results for the correlation 
functions [12, 13]. In the present case, however, the ra- 
pidities can be determined to exponential precision for 
all states such that N s <C N. The remaining N s equa- 
tions for the string centers yield rapidities which in the 
limit L — > oo are quantized as for free particles (full de- 
tails will be published elsewhere [8]). This allows us to 
go analytically much further than in the repulsive case. 

One technical difficulty is that the elements of the norm 
and FF matrices are singular when calculated on exact 
strings with 5 = 0. A more careful analysis however 
shows that all divergences cancel [8] , and that it is possi- 
ble to obtain closed forms for all norms and form factors. 
Single-particle and density FFs between the ground state 
and any other string state have the rather simple forms 

|E*| = ciN\(N -1)1 11 HM/c), (6) 
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| S P| = 0lm(N- 1)1 ]\HM/c), (7) 
where we have defined the fundamental block Hm{x) — 

T ( N-M , - , 2 

These expressions are exact in the limit 



of large L for any finite N [14]. The general expression 
for the norms is 
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Eqs. (6), (7) and (8) give an exact representation of 
each term in the expansion (5) of the desired correlation 
functions for L — > oo for any value of N and c. The 
sums over intermediate states can now be analytically 
performed for various families of excited states, starting 
from the simplest ones. 

Dynamical structure factor. Let us start here with the 
simplest excited states, which are composed of one N- 
string. For arbitrary N, these states give a single coher- 
ent peak, 
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The leading multiparticlc contributions come from two- 
particle N — M : M states, for which we find 

9(^-^ f )fc 4 r 4 (jv) v g^ M (4)^(^f) 
M[ ' J [u,-,S u ]i<*c» h FN - MM ^-^ M ) 

(10) 

where 0(x) is the Heaviside step function, C M — 
2LNi[(N-M)M]i &ndw p M (k) = ^NM(N-M) + ^ is 
the lower threshold for these excitations (there is no up- 
per threshold). p s and /*m are respectively the rapidities 
of the N — M and M strings. Energy and momentum 
conservation enforce pf = jjT [ n ^n-m) 1 * I w _ w m(*0] * > 
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One-body Green function. The most important inter- 
mediate states here are those composed of one N — 1 
string with total momentum fc. For arbitrary N, the 
contribution of these states to the one-body function is 

c*„ , ^/(cL(N-l) 2 ) 

S Q (k,u)= — . (11) 

lla=l ll 1 - 2^J + yac(N-l)) J 

The leading multiparticle contributions come from the 
two-particle TV — 1 — M : M states, 



SZ(k,uj)-- 



e(c-c M )r 4 (iv) 



N-l-M 



(12) 

These form a continuum extending above the boundary 
w M (fc) (obtained from lo p m by replacing N — ► N — 1). 
Energy and momentum conservation constrain p s , pm as 
for S p (replacing TV by N — 1 and ui p M by w M ). 

All other families of excited states can be similarly con- 
sidered, and the desired correlation function can be re- 
constructed to any accuracy required [19]. 

We now consider the special large N limit obtained by 
sending N — > oo while keeping g — cN constant. This 
prescription maintains the binding energy per particle 
Eq/N = —g 2 /12 finite [15]. In this limit, the correlators 
mentioned simplify considerably. For the DSF, we find 
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that the coherent mode contribution becomes 

no,, \ ~ N 2 (irk/a) 2 r/ fc 2 . oN 

SS(k,u) = 2n— 4>J!L-6(u - —). (13) 
sinn TTK/g jv 

and for multiparticle states (the leading M = 1 part), 

TVfc 4 9(^-g 2 /4) 



To certify our results, a number of checks can be made. 
The first is to compute the second density moment (p 2 ), 
which is given by the Hcllmann-Feynman theorem as 
TT ^§i~ ~ — \~ Integrating shows that the one- 
particle part Sq (k, u>) completely saturates this to leading 
order in TV. To go further, we can study the /-sumrule, 
stating that at fixed k the integral ^u)S(k, u) equals 
k 2 N/L. Thus, in the large TV limit, the one TV-string 
states contribute to the /-sumrule as 



fo(k) 
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and therefore to achieve complete saturation, we here 
need higher states (this is natural, since we are computing 
a first frequency moment). The leading M = 1 two- 
particle states contribute to the /-sumrule as 
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and therefore completely saturate the remaining part of 
the /-sumrule when TV string contributions have been 
taken into account. Higher frequency moments further 
suppress the Sq part; in general, the full DSF is therefore 
well approximated by S p (k,oj) = S^k^u) + S p (k, u). 

The static structure factor is the integral of the dynam- 
ical one and can be written as S p (k) — J ^S p (k,u) = 
S p (k) + S p (k), with S p (k) trivial and 
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where ip2(z) is the polygamma function (second deriva- 
tive of the logarithm of the Gamma function). For large 
k the static structure factor is exponentially suppressed 
and is dominated by the two-particle term. The latter 
shows a peak at k/g ~ 1.2 reminiscent of what is obtained 
for the super- Tonks- Girardeau gas- like regime [16]. 
For the one-body function, we similarly find 
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with two-particle states with M = 1 giving the leading 
contribution in the continuum, 
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with F2(x) — 7rx/sinh7rx and x± as defined above. 
To check the relative weight of these states, we use 
the sumrule that the momentum and frequency inte- 
gral of this must equal the density N/L, and obtain 
Coo f iZc t S o(k, w) = f , i-c the one-particle states 
completely saturate this sumrule. However, as we have 
seen with the DSF, calculating higher frequency mo- 
ments would necessitate taking the multiparticle terms 
into account, starting with the M = 1 part. This 
correlation function is therefore well-approximated by 
S 9 (k,w) = S%(k,u) + Sf(k,u). 

The sum rules we have studied are saturated using 
only either the one-particle intermediate states, or the 
simplest two-particle excitations formed by breaking the 
original string in two pieces, one of which is a single ra- 
pidity. For large N, the contributions from higher exci- 
tations are suppressed by progressively higher powers of 
1/N, and can be neglected for most practical purposes. 
The importance of the simplest two-particle excitations 
is highlighted by calculating higher frequency moments. 

Discussion of the results. Both correlation functions 
considered here show similar features. They are domi- 
nated by single-particle coherent contributions, lying on 
lines to — k 2 /N which become dispersionless in the large 
TV limit. There however exist multiparticle continua 
starting from finite mass gap thresholds (u>i ~ g 2 /4 for 
the leading two-particle ones) and extending to arbitrar- 
ily high energies. At the lower threshold, the correlators 
are square-root singular. At high energies, the correla- 
tors decay exponentially, which is reminiscent of massive 
integrable quantum field theories [17]. 

The fact that the flat band single-particle excitations 
account for most of the correlation functions can inter- 
estingly be interpreted as a Mossbauer-like effect, where 
the recoil energy (say upon scattering with a photon in 
Bragg spectroscopy) vanishes because the system acts as 
a single particle of mass TV 3> 1 (the gound state string 
can in fact be viewed as a crystal in rapidity space). It 
should therefore be feasible to observe an analogue of the 
Mossbauer effect on this system. 

The two-particle (N — 1 : 1) part of the DSF is plotted 
in Figure 1 for large N, with g = 1. At the lower thresh- 
old, the DSF diverges as (u — g 2 /4)? . For uj > g 2 /4 the 
DSF is a monotonous decreasing function of u> as long as 
k/g < x c = 1.0565 . . . , whereas for k/g > x c it shows a 
characteristic broad peak, whose position grows like k 2 
for large k and its amplitude decreases like Away 
from this peak the correlator decays exponentially. Four 
fixed momentum cuts are given in Figure 2, showing these 
features in more detail. For the one-body correlator, the 
maximum lies at the lower threshold, with monotonous 
exponential decay at higher energies. 

Concluding remarks. In this letter, we showed 
how intcgrability could be used to compute the zero- 
temperature dynamical structure factor (Eqs (13), (14)) 
and the dynamical one-body correlation function (Eqs 
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FIG. 1: Contribution to the DSF coming from N - 1 : 1 
(two-particle) states, in the large N limit for g = 1. The 
square-root singularity at the lower threshold is accompanied 
by a maximum around fc = [u — gr 2 /4] 2 . 
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FIG. 2: Fixed momentum cuts of Fig. 1, showing the weaken- 
ing of the square-root singularity at higher momentum, and 
the displacement of the maximum towards higher energies. 

(18), (19)) for the ID attractive Bose gas. The cor- 
responding static correlation functions are obtained as 
a subset of our results. The results obtained clearly 
illustrate that this experimentally-realizable quantum- 
mechanical system has many similarities with massive 
integrable quantum field theories, namely the correlation 
functions are given (for fixed k, as a function of to) by 
a coherent peak followed by a multiparticlc continuum. 
Additionally, we have quantified an interesting zero recoil 
energy effect for large particle numbers, which is analo- 
gous to the Mossbauer effect in solids. 

The methods and results obtained in this letter can 
be adapted to describe many other features of the at- 
tractive Lieb-Liniger gas. For example, our expressions 
could lead to a calculation of finite temperature correla- 
tion functions, allowing the description of the quantum 



to classical crossover in this strongly interacting model. 
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